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Abstract 

The main purpose of this study is to develop new thermodynamic models for solution behavior and 
solid-liquid equilibrium in 10 nitrate binary systems of the type 2-1 (Mg(NO,),-H,O, Ca(NO,).- 
H,O, Ba(NO,),-H,O, Sr(NO,),-H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, Al(NO,),-H,O, 
La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,),-H,O) from low to very high concentration at 25 
°C. To construct models, we used different versions of standard molality-based Pitzer approach. To param- 
eterize models, we used all available raw experimental osmotic coefficients data (p) for whole concentra- 
tion range of solutions, and up to supersaturation zone. The predictions of developed models are in excel- 
lent agreement with ¢-data, and with recommendations on activity coefficients (y,) in binary solutions 
from low to very high concentration. The Deliquescence Relative Humidity (DRH), and thermodynamic 
solubility product (as In K°,) of 12 nitrate solid phases, precipitating from saturated binary solutions have 
been calculated. The concentration-independent models for nitrate systems described in this study are of 


high importance for development of strategies and programs for nuclear waste geochemical storage. 
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Introduction 


Computer models that predict solution behavior and solid-liquid-gas equilibria close to 
experimental accuracy have wide applicability. They can simulate the complex changes 
that occur in nature and can replicate conditions that are difficult or expensive to duplicate 
in the laboratory. Such models can be powerful predictive and interpretive tools to study 
the geochemistry of natural waters and mineral deposits, solve environmental problems 
and optimize industrial processes. However, development of comprehensive models for 
natural systems, with their complexity and sensitivity, is a very difficult, time consuming 
and challenging task. The specific interaction approach for describing electrolyte solutions 
to high concentration introduced by Pitzer (1973, 1991) represents a significant advance 
in physical chemistry that has facilitated the construction of accurate computer thermo- 
dynamic models. It was showed that this approach could be expanded to accurately cal- 
culate solubilities in complex brines, and to predict the behavior of natural and industrial 
fluids from very low to very high concentration at standard temperature of 25 °C (Harvie 
et al. 1984; Trendafelov et al. 1995a, 1995b; Christov 1996a, 1998, 1999, 2001, 2002a, 
2002b, 2003a, 2003b, 2003c, 2004, 2005; Christov et al. 1998; Ojkova et al. 1999; Park 
et al. 2009; Kolev et al. 2013; Lach et al. 2018; Guignot et al. 2019; Donchev and Chris- 
tov 2020; Lassin et al. 2020; Donchev et al. 2021; Tsenov et al. 2021), and from 0 to 290 
°C (Petrenko and Pitzer 1997; Christov and Moller 2004ab; Moller et al. 2006, 2007; 
Lassin et al. 2015; Christov 1995, 1996b, 2005, 2007, 2009, 2012, 2020). 

A long term safety assessment of a repository for radioactive waste requires evi- 
dence, that all relevant processes are known and understood, which might have a sig- 
nificant positive or negative impact on its safety. It has to be demonstrated, that the 
initiated chemical reactions don’t lead to an un-due release of radionuclides into the 
environmental geo-, hydro-, and bio-sphere. One key parameter to assess the propaga- 
tion of a radionuclide is its solubility in solutions interacting with the waste. Solubility 
estimations can either be based on experimental data determined at conditions close 
to those in the repository or on thermodynamic calculations. A so called “thermody- 
namic database” created from experimental data is the basis for thermodynamic model 
calculations. Since the disposal of radioactive waste is a task encompassing decades, the 
database is projected to operate on a long-term basis. Chemical models that predict 
equilibrium involving mineral, gas and aqueous phases over a broad range of solution 
compositions and temperatures are useful for studying the interactions between used 
nuclear fuel waste and its surroundings. ‘The reliability of such predictions depends 
largely on the thermodynamic database. Waters of high salinity are not a typical of 
many geochemical environments which may be chosen as future nuclear waste reposi- 
tory sites. This suggests that an accurate description of highly saline waters should be 
required for modeling of chemical interactions in and around nuclear repositories. 
Currently, the most accurate description of saline waters uses the Pitzer ion interac- 
tion model. Extensive thermodynamic databases, which are based on the Pitzer ion 
interaction model was developed within the Yucca Mountain Project (YMTDB: data0. 
ypf.r2) (Sandia National Laboratories 2007), and Thereda project (THermodynamic 
REference DAtabase, THEREDA-Final Report) (Altmaier et al. 2011). Unfortunately, 
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many of introduced in YMTDB and in THEREDA databases Pitzer models are con- 
centration restricted and cannot describe correctly the solid-liquid equilibrium in geo- 
chemical and industrial systems of interest for nuclear waste programs. 

Nitrates are expected to play a significant role in the context of the underground ge- 
ochemical repository of nuclear waste (Lach et al. 2018; Guignot et al. 2019; Donchev 
and Christov 2020; Lassin et al. 2020). More precisely, long-lived, intermediate-level 
radioactive wastes that are planned to be stored in deep clay formations are composed 
of dried sludge from effluent treatments that contain significant quantities of nitrate 
amongst other elements. They are enclosed in specific containers which are placed in 
underground cavities dug in a very low-permeable argillite host rock. “The storage safe- 
ty analyses show that, despite the protection of the concrete or stainless steel-made ex- 
ternal layers of the containers, the formation water of the host rock is likely to migrate 
and reach the waste during the disposal period” (Donchev and Christov 2020; Lassin 
et al. 2020). This would result in the potential dissolution of large amounts of nitrate 
and other elements, resulting in a highly saline, corrosive and oxidative media with a 
high reactivity towards the containment materials and their surroundings, including 
the host rock. Several options for the management of radioactive waste involve their 
preliminary leaching using nitric acid or nitrate salts to recover U and Pu, followed by 
the incorporation of the leaching residues into concrete or metal packages, which are 
then stored underground. For safety analysis purpose, the scenarios envisaged for these 
options assume that formation water returns to the storage compartments sometime 
after the end of the operating period Dossier ANDRA (2005). With a natural pH value 
slightly above 7 in clayey formations, pore water has to flow through basic concrete ma- 
terials before being in contact with the acidic nuclear waste. A large range of pH lead- 
ing to various chemical behaviors can thus be expected in the vicinity of the waste. At 
near-neutral to basic pH values, reactions of hydrolysis, complexation, and formation 
of solid phases can take place and control the fate of radionuclides (Wang et al. 2006; 
Lassin et al. 2020). Therefore, this reactivity must be characterized by development of 
not concentration restricted thermodynamic models, which accurately describe not 
only solution behavior at low molality, but also low and high molality solid-liquid 
phase equilibrium in nitrate systems (Donchev and Christov 2020). The experimental 
data presented in Rard et al. (1977, 2004), Rard and Spedding (1981), Maliutin et al. 
(2020), El Guendouzi and Marouani (2003), and accurate models reported in previ- 
ous studies (Wang et al. 2006; Lach et al. 2018; Guignot et al. 2019; Donchev and 
Christov 2020; Lassin et al. 2020), and in the present work is a step towards this objec- 
tive. It should be noted that THEREDA (Altmaier et al. 2011) do not include models 
for nitrate solutions and solids. The models introduced in YMTDB (Sandia National 
Laboratories 2007), including these for NO,-systems are restricted up to 6 mol.kg". 

In our previous study (Donchev and Christov 2020) we reported very well vali- 
dated accurate thermodynamic models based on Pitzer ion interactions approach for 
7 nitrate binary systems of the type 1-1 (HNO,-H,O, LiNO,-H,O, NaNO,-H,O, 
KNO,-H,O, RbNO,-H,O, CsNO,-H.O, and NH ,NO,-H,O) from low to very high 
concentration at 25 °C. In this study we developed very well validated not concentra- 
tion restricted thermodynamic models for solution behavior and solid-liquid equi- 
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librium in 10 nitrate binary systems of the type 2-1 (Mg(NO,),-H,O, Ca(NO,).- 
H,O, Ba(NO,),-H,O, Sr(NO,),-H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, 
AI(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,),-H,O) from 
low to very high concentration at 25 °C. Models are developed on the basis of Pitzer 
ion interactions approach. The models for nitrate systems described in this study are 
of high importance, especially in development of strategies and programs for nuclear 
waste geochemical storage. These models are also of interest for industrial application, 
such as production and purification of nitrate compounds. 


Methodology 


The models for nitrate binary systems have been developed on the basis of Pitzer’s semi- 
empirical equations (Pitzer 1973, 1991). Since the Pitzer’s representation of the aque- 
ous phase is based on the excess free energy, all the activity expressions are consistent, 
allowing different kinds of data (e.g., osmotic, emf, and solubility measurements) to 
be used in the parameter evaluations and other thermodynamic functions to be calcu- 
lated (Christov and Moller 2004a, Christov and Moller 2004b; Christov 2007, 2009, 
2012). Pitzer approach has found extensive use in the modeling of the thermodynamic 
properties of aqueous electrolyte solutions. Several extensive parameter databases have 
been reported. These include: 25 °C database of Pitzer and Mayorga (1973, 1974) 
(summarized also in Pitzer 1991); of Kim and Frederick (1988); YMTDB (Sandia 
National Laboratories 2007), and THEREDA (Altmaier et al. 2011). However, some 
of the models in all of these databases are concentration restricted, and do not include 
all minerals precipitating from saturated and supersaturated binary and mixed systems. 
The most widely used are databases of Chemical Modelling Group at UCSD (Uni- 
versity California San Diego): at 25 °C (Harvie et al. 1984; Park et al. 2009), and T- 
variation (from 0 to 300 °C) (Christov and Moller 2004a, Christov and Moller 2004b; 
Moller et al. 2006, 2007; Christov 2009). Some of comprehensive minerals solubilitiy 
YMTDB (Sandia National Laboratories 2007), and THEREDA (Altmaier et al. 2011) 
databases also contain concentration restricted models for some low-, or high- concen- 
tration binary and mixed sub-systems with strong association reactions in unsaturated 
solutions. The concentration restricted sub-models are developed using experimental 
activity data in binary solutions, and solubility data in binary and high order systems 
up to maximum concentration (m(max)), which is much lower than concentration 
of saturated or supersaturated binary and mixed solutions (m(sat)). Such a restricted 
models predict minerals solubility, which is in pure agreement with experimental data. 


The Pitzer’s equations 


According to Pitzer theory electrolytes are completely dissociated and in the solution 
there are only ions interacting with one to another (Pitzer 1973; Pitzer and Mayorga 
1973). Two kinds of interactions are observed: (i) specific Coulomb interaction be- 
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tween distant ions of different signs, and (ii) nonspecific short-range interaction be- 
tween two and three ions. The first kind of interaction is described by an equation of 
the type of the Debye-Hueckel equations. Short-range interactions in a binary system 
(MX(aq)) are determined by Pitzer using the binary parameters of ionic interactions 
(8,8, C*, and B®). The Pitzer’s equations (1 to 4) are described and widely discussed 
in the literature (Harvie et al. 1984; Moller et al. 2006, 2007; Christov and Moller 
2004a, Christov and Moller 2004b; Christov 2005). Here only the expression for the 
activity coefficient of the interaction of cation (M) with other solutes, 


Your is given: 


Inv, =2,,F +> m, (2By,(1)+ZCy,)+ > m, [29 LIne | LE MI + [Zu] >) >) m.11,Cog 


ve m, (24a )+ 2, Dial nal (1) 


Equation (1) is symmetric for anions. The subscripts c and a in eqn 1 refer to 
cations and anions, and m is their molality; z is the charge of the M* ion. B and ® 
represent measurable combinations of the second virial coefficients; C and ) represent 
measurable combinations of third virial coefficients. B and C are parameterized from 
single electrolyte data, and @ and ¢) are parameterized from mixed solution data. The 
function F is the sum of the Debye-Hueckel term, 


-A? [VI/ (1+ bV I) + (2/b) (n(1 + bv I) ], (2) 


and terms with the derivatives of the second virial coefficients with respect to ionic 
strength (see Harvie et al. 1984). In Eq. (2), b is a universal empirical constant assigned 
to be equal to 1.2. A? (Debye-Hiickel limiting law slope for the osmotic coefficient) is 
a function of temperature, density and the dielectric constant of water (Christov and 
Moller 2004b). 

For the interaction of any cation M and any anion X in a binary system MX-H,O, 
Pitzer assumes that in Eq. (1) B has the ionic strength dependent form: 


Bee = BOL er Bie o(a, VI) + (3) 

+8,  g(a,VD, (3A) 

where g(x) = 2[1 - (1 + x)e*] / x’ with x = aV lor a,V I. x terms are function of elec- 
trolyte type and does not vary with concentration or temperature. 


In Eq. 1, the ® terms account for interactions between two ions i and j of like 
charges. In the expression for ®, 


®, = 0, + =O (1), (4) 


0,, is the only adjustable parameter. ‘The "0. (I) term accounts for electrostatic unsym- 
metric mixing effects that depend only on the charges of ions i and j and the total ionic 
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strength. ‘The ¢,,, parameters are used for each triple ion interaction where the ions are 
not all of the same sign. Their inclusion is generally important for describing solubili- 
ties in concentrated multicomponent systems. Therefore, according to the basic Pitzer 
equations, at constant temperature and pressure, the solution model parameters to be 
evaluated are: 1) pure electrolyte 8, 8“, and C® for each cation-anion pair; 2) mixing 
0 for each unlike cation-cation or anion-anion pair; 3) mixing ) for each triple ion 
interaction where the ions are all not of the same sign. 

Fluids commonly encountered in natural systems include dissolved neutral species 
(such as carbon dioxide (CO; 3 SiO 544) and Al(OH),°(aq)). To account neutral specie 
interactions in aqueous solutions the UCSD Chemical Modelling Group included in 
their models additional terms to Pitzer equations, denoted as ese or hy Re and ¢ 


(Eq. (1)) (Harvie et al. 1984; Moller et al. 2006, 2007). 


N, A,X. 


The 8” parameter (Eqn. 3A) for 2—2 type of electrolytes 


Pitzer and Mayorga (1973) did not present analysis for any 2—2 (e.g. MgSO,-H,O) or 
higher {e.g. 3-2: Al,(SO,),-H,O} electrolytes. Indeed, they found that three 6, 6, 
and C? parameters approach (see Eqns. 1 and 3) could not accurately fit the activity data 
for these types of solutions. For these electrolytes mean activity (y,) and osmotic (y) 
coefficients drop very sharply in dilute solutions, while showing a very gradual increase, 
with a very wide minimum at intermediate concentration. Pitzer concluded that this 
behaviour is due to ion association reactions and that the standard approach with three 
evaluated solution parameters cannot reproduce this behaviour. This lead to a further 
(Pitzer and Mayorga 1974) modification to the original equations for the description of 
binary solutions: parameter 8*(M,X), and an associated a,V I term are added to the B,,, 
expression (see Eqn. (3A). Pitzer presented these parameterizations assuming that the 
form of the functions (i.e. 3 or 4 8 and C ® values, as well as the values of the « terms) 
vary with electrolyte type. For binary electrolyte solutions in which either the cationic 
or anionic species are univalent (e.g. NaCl, Na,SO,, or MgCl), the standard Pitzer 
approach use 3 parameters (i.e. omit the 8 term) and «, is equal to 2.0. For 2—2 type 
of electrolytes the model includes the 8 parameter and «, equals to 1.4 and «, equals 
to 12. This approach provides accurate models for many 2—2 binary sulfate (Pitzer and 
Mayorga 1974; Christov 1999, 2003a) and selenate (Christov 2003a; Christov et al. 
1998) electrolytes, giving excellent representation of activity data covering the entire 
concentration range from low molality up to saturation and beyond. 


Inclusion of “standard Pitzer approach” 8” parameter into a models for I-I, 
2-1, 3-1, 4-1, 1-2, 1-3, and 3—2 type of electrolytes 


Some authors found that there are some restrictions limited the potential of the model 
to describe correctly activity and solubility properties in some binary electrolyte sys- 
tems with minimum one univalent ion (see Petrenko and Pitzer 1997, 2012; Gruszk- 
iewicz and Simonson 2005; Lach et al. 2018; Guignot et al. 2019; Lassin et al. 2020), 
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and of 3—2 type (see Christov, 2002ab, 2003b) at very high molality using classical 3 
parameters (8, 8", and C*) approach. According to discussion in Christov (2004, 
2005, 2012) and in Lassin et al. (2015), there is one major factor which determined 
these restrictions: type of ¢ (osmotic coefficient) vs. m, or y, (activity coefficient) vs. 
m dependences at high concentration. For all these systems y vs. m, or y, vs. m curves 
have a wide maximum at molality approaching molality of saturation: “LiCl(aq) type”: 
see Lassin et al. (2015), “FeCl,(aq) and FeCl,(aq) type”: see Christov (2003c, 2004); 
“HNO, (aq) type”: see Donchev and Christov (2020);“Al (SO,),(aq), Cr,(SO,),(aq) 
type”: see Christov (2002a, 2002b, 2003b). 

To describe the high concentration solution behaviour of systems showing a 
“smooth” maximum on y, vs. m dependence, and to account strong association reac- 
tions at high molality, Christov (1996a, 1998, 1999, 2001, 2005) used a very simple 
modelling technology: introducing into a model a fourth ion interaction parameter 
from basic Pitzer theory {8 in Eqn. (3A)}, and varying the values of «, and «, terms 
in Eqs. (3 and 3A))). The author also found that by variation of the values of «, and 
a, terms it is possible to vary the concentration range of binary solutions at which as- 
sociation reactions become more important and should be account by introducing 8” 
parameter. According to Christov (2005), model which uses «, = 1.4 and a, = 12 ac- 
counts association only at low molality solutions (see also Christov and Moller (2004b) 
for Ca(OH),-H,O model). According to previous studies of one of the authors (Chris- 
tov) an approach with 4 ion interaction parameters (8,8, 8®,and C*), and accept- 
ing a, = 2, and varying in , values can be used for solutions for which ion association 
occurs in high molality region. This approach was used for binary electrolyte systems 
of different type: 1-1 type {such as HNO,-H,O, LiNO,-H,O (Donchev and Chris- 
tov 2020), CsF -H,O (Tsenov et al. 2021), and LiCl-H,O (Lassin et al. 2015)}, 2-1 
{such as NiCL-H,O, CuCl,-H,O, MnCl -H,0O, CoCL,-H,O: (Christov 199G6a, 1999); 
FeCl,-H,O: (Christov 2004); Ca(NO,),-H,O: (Lach et al. 2018); UO,(NO,),-H20 
(Lassin et al. 2020)}, 1-2 {such as Na,Cr,O,-H,O: (Christov 2001); K,Cr,O_-H,O: 
(Christov 1998)}, 3-1 {such as FeCl,-H,O: (Christov 2004); Ln(NO,),(aq): (Guignot 
et al. 2019)}, and 3-2 {such as Al (SO,),-H,O, Cr,(SO,),-H,O, and Fe,(SO,),-H,O: 
(Christov 2001, 2002a, 2003b, 2004, 2005)}. The resulting models reduce the sigma 
values of fit of experimental activity data, and extend the application range of models 
for binary systems to the highest molality, close or equal to molality of saturation 
{m(sat)}, and in case of data availability: up to supersaturation. For example, aque- 
ous complexes free 4 parameters model for LiCl-H,O system predicts LiCl.nH2O(s) 
solubilities from 0 to 200°C and up to 40 mol.kg" (Lassin et al. 2015). The resulting 
accurate 4 - parameters solution models are used directly to determine InK°, values of 
precipitated solid phases using solubility approach (Harvie et al. 1984; Christov 1995, 
Christov 1996a, Christov 1996b, Christov 2005, Christov 2012; Christov and Moller 
2004a, Christov and Moller 2004b). Therefore, the developed not high molality re- 
stricted parameterization, were used without any changes for development of solid- 
liquid equilibrium models for high order systems. ‘Thus, models for Al,(SO,,),(aq) and 


Cr,(SO,),(aq) are used without additional adjustments to construct a model for multi- 


396 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) 


component (Na+K+NH,+Mg+Al+Cr+SO ,+H,O) system (Christov, 2002ab, 2003b). 
Four parameters (8, 6, 8 and C*) models for NiCl,(aq), CuCl,(aq), MnCl (aq), 
and CoCl,(aq) are used for construction of Na-K-Rb-Cs-Ni-Co-Cu-Mn-Cl-H,O 
model (Christov 1996a, 1999). Four parameters models for FeCl,(aq)and FeCl, (aq) 
are directly used in development of high accuracy minerals solubility model for (Na 
+K+Mg+Fe(II)+Fe(II)+Cl+SO,+H,O) system (Christov 2004). A model for binary 
systems Na,Cr,O,(aq), K,Cr,O,(aq) was used without any changes to develop a com- 
prehensive model for: (Nat+K+Cl+SO,+Cr,O,+H,O) system (Christov 1998, 2001), 
and Ca(OH),(aq) model is used as s strong base for H+Na+K+Ca+OH+Cl+SO ,+H,O 
model from from 0 to 300°C (Christov and Moller, 2004b). 


Results and discussions 


Model parameterization and validation of models for binary 2—1, 3—I,and 4—| 
type nitrate systems 


In this study we developed new thermodynamic models for solution behavior and 
solid-liquid equilibrium in 10 nitrate binary systems of the type 2-1 (Mg(NO,).- 
H,O, Ca(NO,),-H,O, Ba(NO,),-H,O, Sr(NO,).-H,O, and UO,(NO,),-H,O), 
3-1 (Cr(NO,),-H,O, AI(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 
(Th(NO,) ,-H,O) from low to very high concentration at 298.15 K. New sets of Pitzer 
ion interaction binary parameters are evaluated using available raw experimental osmot- 
ic coefhicients (p) data for whole molality range of solutions. Rard and co-authors (1977, 
1981) reported an extensive experimental activity database for rare earth nitrate systems. 
Data of Rard and Spedding (1981) are used to parametrize the model for Lu(NO,),- 
HO system. The 9 vs. m data for remaining 9 nitrate solutions under study are given in 
Mikulin (1968), and Robinson and Stokes (1959). Reference » vs. m data sets of Miku- 
lin (1968), and Robinson and Stokes (1959) are in a good agreement. Data of Mikulin 
(1968) for La(NO,),-H,O are also in good agreement with the data of Rard (1987). 
However, the data of Mikulin (1968) cover the whole molality range of unsaturated 
and saturated solutions. In case of Ca(NO,),-H,O, UO,(NO,),-H,O, and Th(NO,),- 
HO systems, Mikulin also reported data for supersaturated solutions. In this study we 
parameterize the models using 1) all data of Mikulin (1968) for whole molalty range of 
unsaturated solutions from 0.1 m to m(max), 2) the data points at saturation (¢(sat)) 
(from Mikulin 1968), and 3) data for supersaturated Ca(NO,),-H,O, UO,(NO,).- 
HO, and Th(NO,),-H,O solutions (from Mikulin 1968), and 4) all experimental and 
recommended data of Rard and Spedding (1981) for Lu(NO,),-H,O system. 

In parameterization we used the value of Debye-Hiickel term (A®) equals to 0.39147 
(Christov 2007, 2009, 2012). Following the parameterization scheme described in pre- 
vious paragraph the model for all 10 binary nitrate solutions is parameterized using 
two different approaches: (I) standard for N-1 electrolytes (N = 2, 3, or 4) approach 


with 3 ion interaction binary parameters (8, 6, and C’) and setting «, term equals 
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to 2, and «, = 0.0, and (II) an extended approach with four Pitzer ion interaction bi- 
nary parameters (8, 8", 8°,and C*) and varying in the values of «, and a, terms. As 
a first step in parameterization we used classical 3 parameters approach (I) and evaluate 
binary parameters using all available raw » data for whole molality range of solutions. 
As a next step, using the same ¢ data we re-parameterize the models on the basis of 
extended approach (II), and using three « — combinations: (Ia) «, = 2 and «, = 1, and 
(IIb) «, = 2 and «, = -1 (Christov 1996a, 1998, 1999, 2004, 2005) and (IIc) «, = 2 and 
a, = 0.3 (Guignot et al. 2019; Donchev and Christov 2020; Donchev et al. 2021). It 
was found that more combinations in “alfa” values do not improve the fit of data used 
in parameterization. The main criterion in the choice of established parameterization 
was the value of standard deviation (co) of fit of used ¢ data, i.e. parameterization with 
the lowest sigma value is accepted. For definition of sigma (co) see Christov and Moller 
(2004b), and Christov (2007, 2009, 2012). It was found that for 2 of studied systems 
Ba(NO,),-H,O, and Sr(NO,),-H,O, the approach (I) with 3 parameters (8, 6, C) 
give an acceptable agreement with the data. For these systems introducing into a model 
of fourth (8) parameter do not improve considerably the fit of data. For all other ni- 
trate systems under study we construct a model on the basis of extended approach (I), 
and using different combinations of “alfa’values: for Lu(NO,),-H,O system. 1) «, = 2 
and a, = 0.3 (approach IIc), and 2) «, = 2 and a, = 1 (approach IIa), and 3) «, = 2 and 
a, = -1 (approach Ib). The resulting models fits the data up to supersaturation zone 
(m(max) = 14.77 m in Ca(NO,),-H,O) with sigma values, which is much less than the 
sigma values of models of Pitzer and Mayorga (1973), and of Kim and Frederick (1988). 

On next Figure 1 we present a comparison of osmotic coefhcients in nitrate bi- 
nary solutions 2-1 (Mg(NO,),-H,O, Ca(NO,),-H,O, Ba(NO,),-H,O, Sr(NO,).- 
H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, Al(NO,),-H,O, La(NO,),-H,O, 
Lu(NO,),-H,O), and 4-1 (Th(NO,) -H,O) calculated by the accepted models devel- 
oped here (heavy solid lines), and models developed by other authors (dashed lines 
and light solid lines: Pitzer and Mayorga (1973), Kim and Frederick (1988), Rard and 
Spedding (1981) (for Lu(NO,),-H,O system only), Rard et al. (2004) (for Mg(NO,).- 
H,O system only), and Wijesinghe and Rard (2005) (for Ca(NO,),-H,O system only). 
The recommended osmotic coefficients values given in literature at 25 °C are given on 
Fig. 1 by symbols. The vertical lines on the figures denote the molality of solutions 
saturated with corresponding nitrate solid phase (m(sat)), taken from Mikulin (1968). 
Excellent accepted new model (heavy solid line) — experiment (symbols) agreement 
has been obtained for all 10 systems and from low (see figures for Mg(NO,),-H,O, 
Ca(NO,),-H,O systems) and up to very high molality. As is shown Fig. 1, the new 
model for Ca(NO,),-H,O is in excellent agreement not only with the data at high mo- 
lality ((m(max) = 14.77 m), but contrary to the models of Kim and Frederick (1988) 
and Wijesinghe and Rard (2005) also in low molality range. It should be noted that 
all reference models presented on Fig. 1 by dashed, dashed-dotted, and light solid lines 
(Kim and Frederick (1988), Pitzer and Mayorga (1973), Rard et al. (2004), Wijesinghe 
and Rard (2005), and Rard and Spedding (1981)) have been constructed on the basis 


of standard Pitzer approach with 3 interaction parameters. Therefore, these models 
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Figure |. Comparison of model calculated (lines) osmotic coefficients (y) of Mg(NO,), Ca(NO,),, 
Ba(NO,),, Sr(NO,), UO,(NO,),, Cr(NO,) i AI(NO,) * La(NO,),, Lu(NO,),, and Th(NO,), in binary so- 
lutions 2-1 (Mg(NO,).-H,O, Ca(NO,),-H,O, Ba(NO,),-H,O, Sr(NO,).-H,O, and UO,(NO,),-H,O), 
3-1 (Cr(NO,),-H,O, AI(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,),-H,O) 
against molality at T = 298.15 K, with recommendations in literature (symbols). For Mg(NO,),-H,O 
and Ca(NO,),-H,O systems an enlargement of the low molality corner is also given. Heavy solid lines 
represent the predictions of the developed in this study and accepted models. Dashed-dotted, dashed 
and light solid lines represent the predictions of the reference models of Kim and Frederick (1988), of 
Pitzer and Mayorga (1973), of Rard et al. 2004 (for Mg(NO,),.-H,O), of Wijesinghe and Rard (2005) 
(for Ca(NO,),-H,O), and of Rard and Spedding (1981) (for Lu(NO,),-H,O). On Figures YM denotes 
YMTDB (Sandia National Laboratories 2007). For Lu(NO,),-H,O the experimental data and recom- 
mended data are taken from Rard et al. (1977) (open squares), and Rard and Spedding (1981) (crosses), 
respectively. For all other systems the experimental data of Mikulin (1968) are used (open squares and 
open triangles). The molality of stable and metastable (for Ca(NO,),-H,O) crystallization of solid nitrate 
phases (m(sat)) is given on all figures by vertical lines (see Table 1). 
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Figure |. Continue. 


cannot reproduce well the experimental data (Fig. 1). To illustrate this conclusion on 
Fig. 1 we give the predictions of two new models for La(NO,),-H,O system. As it is 
shown the 3 parameters model (light solid line) is in pure agreement with the data. 

The models for all nitrate binary systems under study are also validated by com- 
parison with recommendations given in literature (Rard and Spedding (1981) for 
Lu(NO,),-H,O); and Mikulin (1968)) (for all other 9 systems under study) on the 
mean activity coefficients (y,). These recommendations on y, are model-dependent. 
Therefore, they are not used in parameterization process, and only to validate the re- 
sulting models. The comparisons between predictions of new developed models and 
reference recommendations, which are not given here, show an excellent agreement 
from low to very high concentrations. 


Deliquescence relative humidity (DRH) calculations 


Deliquescence of single inorganic salt or their mixture is a process of spontaneous 
solid-liquid phase change. It is a process in which a soluble solid substance sorbs water 
vapor from the air to form a thermodynamically stable saturated aqueous solution on 
the surface of the particle. It is occurring when relative humidity (RH) in the gas-phase 
environment is at, or above deliquescence relative humidity (DRH) of the salt, or 
mutual deliquescence relative humidity (MDRH) of a salt mixture. Within the solid- 
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liquid equilibrium model, relative humidity is related to water activity(z,) (Clegg et 
al. 1998; Christov 2009, 2012; Donchev and Christov 2020) according to Eqn. (5): 


a,=P./P° = RH/100, (5) 


where P, and P° are the vapor pressure of the saturation solution and pure water, re- 
spectively, at given temperature. As a result, both DRH and MDRH of saturated surface 
solutions depend of temperature, the salt stoichiometry, and the solution composition. 
This process is of interest in many areas, such as heterogeneous chemistry of inorganic 
salts, corrosion of metals in wet atmosphere, in studies of chemistry of sea-type aero- 
sol atmospheric system (Kolev et al. 2013), and especially in development of strategies 
and programs for nuclear waste geochemical storage. Because of very high complicity 
of experiments, the relative humidity DRH experimental data are sparse. ‘Therefore, 
different sophisticated thermodynamic models have been proposed and developed to 
describe the deliquescence behavior of inorganic salts at wet conditions. In our previ- 
ous studies it was showed that calculations based on not high concentration restricted 
Pitzer models can be used for accurate determinations of both DRH and MDRH of 
saturated solutions in a wide range of temperatures, and compositions (Christov 2009, 
2012; Donchev and Christov 2020). On the basis of evaluated binary parameters (8, 
8, B,and C*) in this study we also determine water activity (4,) and Deliquescence 
Relative Humidity (DRH (%)) (eqn. 5) of 12 solid phases crystallizing from saturated 
binary nitrate solutions [Mg(NO,),.6H,O(s), Ca(NO,),.4H,O(s), Ca(NO,),.3H,O\(s), 
Ba(NO,),(s), Sr(NO,),(s), UO,(NO,),.6H,O(s), AI(NO,),.9H,O(s), Cr(NO,),(s), 
La(NO,),.6H,O(s), La(NO,),(s) Lu(NO,),.5H,O(s) and Th(NO,),.6H,O(s)]. Note that 
the widely used databases of Pitzer (1991), Pitzer and Mayorga (1973), Pitzer and Kim 
(1974) and Kim and Frederick (1988) do not consider solid phases. The results of calcu- 
lations are given in Table 1. The model DRH predictions are in excellent agreement with 
the experimental data determined using isopiestic method, and given in Mikulin (1968). 
According to model calculations the solid-liquid phase change of Ca(NO,),.3H,O(s), 
and Lu(NO,),.5H,O(s) occurs at lowest relative humidity of environment. It can be con- 
cluded that the solid-liquid phase change of solid nitrates of Lanthanide metals is more 
activated in the presence of calcium in the nuclear storage environment. 


Determination of thermodynamic solubility product (K°.,) of precipitates 


In this study we determine the thermodynamic solubility products (as K° ) of solid 
phases, precipitating from saturated nitrate binary solutions, s.a. anhydrous Ba(NO,),(s) 
and hydrate Ca(NO,),.3H,O(s), precipitating in Ba(NO,),-H,O and Ca(NO,),-H,O. 
The K°, have been determined on the basis of evaluated binary parameters and using 
experimental m(sat) solubility data, and using the following relationships (Christov 
2005, 2007, 2009, 2012): 


K°sp (Ba(NO,),) = 4. Y (sat)? . m(sat)? 
K°sp(Ca(NO,),.3H,O) = 4. Y (sat) >, m(sat) °. a (sat) * (G) 
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Table |. Comparison between model calculated and recommended values of the Deliquescence Relative 
Humidity [DRH (%) = a, (sat). 100; where a, (sat) is activity of water at saturation] and of the logarithm 
of the thermodynamic solubility product (as InK”_) of nitrate solid phases crystallizing from saturated 
binary solutions at T = 25° C. 


Salt composition m(sat) (exp) Ink iff DRH(%) 
(mol.kg’) This work Reference data This work Reference 
calculated calculated data* 
Mg(NO,),.6H,O(cr) 5.06 7.0098 7.02% 52.32 52.90 
Ca(NO,) , 4H,O(cr) (stable solid) 8.41? 4.4362 4.53 49.07 49.10 
Ca(NO,) , 3H,O(cr) (metastable solid) 14.77? 6.6449 5.34° (m/(sat) = 15.0 m) 22.52 - 
Ba(NO,),(cr) 0.39? -5.125 - 98.61 98.60 
St(NO,), .4H,O (cr) 3.76% 0.0327 - 84.83 84.80 
UO,(NO,), .6H,O(cr) 3.21 5.3022 5.251° 73.44 73.60 
AI(NO,),.9H,O (cr) 3.16 4.3081 : 59.88 60.20 
Cr(NO,), (cr) ¢ 1.4° 1.2097 - 86.38 - 
La(NO,), .6H,O (cr)¢ 4.6154 2.1599 297° 62.38 - 
La(NO,), (cr) * 2.94# 1.4704 < 77.73 77.60 
Lu(NO,),5H,O (cr) 6.8154 10.7681 10.67 31.27 - 
Th(NO,),.6H,O(cr) 4.00 4.4886 4,71° (as Th(NO,),.5H,O) 54.46 55.0 


Guignot et al. (2019); ‘Accepted m(sat) molality and stoichiometry of solid phase. 

As a next step, using the accepted new developed parameterizations, and experi- 
mentally determined molalities (m(sat) of the saturated binary solutions (Mikulin 
1968; Guignot et al. 2019; Lassin et al. 2020)) we calculate the logarithm of the 
thermodynamic solubility product (In K°) of twelve nitrate solid phases crystal- 
lizing from saturated binary nitrate solutions at 25 °C (Eqn. (6)). The model cal- 
culations are given in Table 1. With only 2 exceptions (for Ca(NO,),.3H,O(s) and 
La(NO,),.6H,O(s)) a good agreement has been obtained with calculations of Lach 
et al. (2018), Lassin et al. (2020), and Guignot et al. (2019) for all nitrate solids. The 
In K°. | differences are mainly due on the 1) different m(sat) values used in calcula- 
tions (see Egn. (6)), and 2) different experimental data source with different m(max) 
values used in parameterization. 


Summary and conclusions 


In this study we developed new thermodynamic models for solution behavior and 
solid-liquid equilibrium in 10 nitrate binary systems of the type 2-1 (Mg(NO,).- 
H,O, Ca(NO,),-H,O, Ba(NO,),-H,O, Sr(NO,),-H,O, and UO,(NO,),-H,O), 
3-1 (Cr(NO,),-H,O, Al(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 
(Th(NO,),-H,O) from low to very high concentration at 25 °C. To parameterize models 
for binary systems we used all available raw experimental osmotic coefficients data () 
for whole concentration range of solutions, and up to saturation point. Data for super- 
saturation zone, available for Ca(NO,),-H,O, UO,(NO,),-H,O, and Th(NO,),-H,O 
systems, are also included in parameterization. To construct models, we used different 
versions of standard molality-based Pitzer approach. It was established that with only 2 
exceptions (Ba(NO,),-H,O, and UO,(NO,),-H,O) application of extended approach 
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with 4 parameters (8, 6", 8°,and C*) and variation of «, term in fundamental Pitzer 
equations leads to the lowest values of standard model-experiment deviation. The pre- 
dictions of new developed here models are in excellent agreement with experimental os- 
motic coeflicients data (see Fig. 1), and with recommendations on activity coefficients 
(not given here) in binary solutions from low to very high concentration: up to 14.77 mol. 
kg" in Ca(NO,),-H,O. The Deliquescence Relative Humidity (DRH), and thermody- 
namic salubiliey reduce (as In K°, ») of of 12 solid phases crystallizing from saturated bi- 
nary nitrate solutions [Mg(NO, ), Hae sO(s), Ca(NO,),.4H,O(s), Ca(NO,),.3H,O\(s), 

Ba(NO,),(s), Sr(NO,),(s), UO (NO, ),-6H,O(s), AL(NO, Jeet COKS) Cx(NO, ),(s), 

La(NO,),.6H,O(s), La(NO,),(s) Lu(NO,),.5H,0() and “Th(NO,),.6H,O(s)] fae 
been determined on the basis of evaluated binary parameters and using experimental 
m(sat) solubility data. Model predictions are in good agreement with available refer- 
ence data. The accurate solid-liquid equilibrium models for nitrate systems described 
in this study are of high importance for development of strategies and programs for 
nuclear waste geochemical storage. 
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